RANGE ADAPTABLE SYSTEM FOR DETERMINING THE 
ANGULAR POSITION AND DISTANCE OF A RADIATING 
POINT SOURCE AND METHOD OF EMPLOYING 

BACKGROUND OF THE INVENTION 

The present invention relates to methods of employing a range 
adaptable system for determining the angular position and distance 
of a radiating point source as well as the system itself. In the 
present invention, radiation from a point source travels through a 
window mask pattern and onto the surface of a sensor. In the prior 
art, this general concept is well known, however, the prior art 
fails to teach or suggest the specific techniques employed by 
Applicant to accurately measure angular position and distance using 
such a structure. 

Prior art systems have used masked or coded apertures placed 
above multi-element detectors. All of these techniques have 
attempted to determine the position of a projected pattern on a 
detector surface by comparing signals from fixed detector elements 
arranged in a unique pattern or from comparison to fixed reference 
signals previously stored in a computer. Such techniques result in 
limited resolution due to the fact that the magnified mask 
projection due to point source illumination does not precisely 
correlate or match the fixed detector pattern or reference signal. 

As an example of the effect of magnification on the alignment 
process, resolving patterns with as many as 100 line/space pairs, 
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mask and image is used to determine a lag function of the received 
mask pattern. A mosaic pattern is disclosed as yielding results 
superior to those yielded through the use of a random array. In 
the Fenimore x 633 patent, the mask and detector are either one or 
two dimensional. The Fenimore '780 patent only discloses two 
dimensional imaging. The present invention differs from the 
teachings of the Fenimore patents since it does not compare image 
to a fixed reference pattern. 

U.S. Patent No. 4,435,838 to Gourlay discloses the concept of 
imaging planes of various distances in the field of tomographic 
imaging. The Gourlay technique improves the prior art process of 
selecting various image sizes to correlate with detector response 
to select a particular depth plane. Gourlay teaches the setting of 
a coded aperture at a distance D/d, where D is the object-to- 
detector distance and d is the object-to-mask distance, so that the 
image can be kept at the same magnification, thereby simplifying 
the correlation process. The present invention differs from both 
that of Gourlay and that of prior art systems in the field of 
tomographic imaging, since, in the present invention, various 
object-to-mask distances or various image sizes are not chosen to 
select a particular depth, but, instead, the mask image scale is 
automatically determined using frequency domain techniques. 

Applicant is also aware of U.S. Patent Nos. 5,408,323 to 
Mitchelson and 6,141,104 to Schulz. Each of these references 
teaches the concept of the use of a random pattern mask or 



there must be less than one part in 400 magnification in order to 
preserve correlation with a fixed reference. This level of 
magnification may typically occur at a distance of between 2 and 4 
meters (6 to 13 feet), which represents the outer range of most 
optical position measurement systems. In order to facilitate 
operation at closer distances, with high resolution performance, a 
means of scale identification is required. It is with this notion 
in mind, among others, that the present invention was developed. 
The following prior art is known to Applicant: 
U.S. Patent Nos. 4,193,689 and 4,314,761, both to Reymond et 
al. disclose arrangements for locating radiating point sources 
including the use of a single axis detector array and an aperture 
mask containing a slit to project light onto the array. In the 
earlier patent, a cylinder lens is used, and in the later patent, 
the cylinder lens is replaced with an aperture that can be 
shuttered to allow light in from preferred directions. In each 
case, three single axis arrays are used, and three cameras are 
required to compute three planes that intersect to define a point 
in space. The embodiments of the present invention differ from the 
teachings of the Reymond et al. patents since they include 
computing means to determine scales and shifts of image components 
and wherein the mask is not a single slit. 

U.S. Patent No. 4,209,780 to Fenimore et al. and U.S. Patent 
No. 4,389,633 to Fenimore disclose the technique of using uniformly 
redundant arrays for coded aperture imaging. Correlation between 



universal redundant array. In each of these patents, measurement 
is conducted by comparing actual measurements to previously stored 
reference measurements. Over and above the reasons why the present 
invention, in its embodiments, patentably distinguishes from the 
Reymond et al. patents, the present invention distinguishes from 
Fenimore, Mitchelson, Gourlay and Schulz as specifically avoiding 
comparison of actual measurements to reference measurements. 

SUMMARY OF THE INVENTION 
The present invention relates to methods of employing a range 
adaptable system for determining angular position and distance 
of a radiating point source as well as the system itself. The 
present invention includes the following interrelated objects, 
aspects and features: 

(1) The present invention is disclosed in terms of two 
preferred embodiments. In each of these embodiments, a patterned 
mask or coded aperture is located at a distance from a linear array 
comprising a detector surface. The mask pattern consists of a 
transmissivity pattern varying in one dimension, e.g., in the Y- 
axis. In each of the embodiments, a point radiating source 
illuminates the aperture to cast an image onto the array. 
Computing means is employed to identify one or more frequencies in 
the frequency domain to determine the image scale and shift along 
the detector array axis. 
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(2) In further amplification of the BACKGROUND OF THE 
INVENTION/ the main distinction between the embodiments of the 
present invention and the prior art is that, in the present 
invention, a method of determination is employed that automatically 
allows determination of the magnification of the aperture image, 
with this method employing frequency domain techniques. The prior 
art fails to teach or suggest a way to account for image 
magnification that occurs due to the proximity of the source and 
the detector. Through the use of the techniques disclosed herein, 
the aperture pattern is automatically re-scaled to match that of 
the actual image, so that an accurate determination of pattern 
shift can be made. 

(3) A first embodiment of the present invention has two 
variations, one of which employs the use of multiple single 
frequency components and phase methodology, the second of which 
uses multiple single frequency components as well as a variable 
frequency component, all of which are identified in "frequency 
domain." The particular locations of frequency components as well 
as the corresponding phases of single frequency components allow 
the determination of the magnified image shift on the detector 
axis. The flowchart of Figure 16 depicts these variations. 

(4) In the first variation, measurement of the frequency and 
phase of the lowest frequency component is used to determine a 
coarse position measurement and is also used to measure the 
particular cycle of the next higher frequency component. In the 



same manner, the frequency and phase of the next higher frequency 
component is used to determine the detector position to a finer 
degree. As before, at the same time, the next higher frequency 
component is also used to determine the particular cycle of the 
next higher frequency component as the measurements become finer 
and more accurate consecutively. 

(5) Eventually, the exact pattern shift is determined and from 
this information, the angle of incidence between the point source 
of radiation and the detector is determined. 

(6) One example of an application uses a 4, 096 element CCD 
array with a mask pattern containing four frequencies, including 
the fundamental frequency f 0 with period of 4,096 pixels, ^=4*^, 
f 2 =16*f 0 , and f 3 =64*f 0 . The mask length (-8192 pixels) is longer 
than the detector to ensure that the detector is fully covered by 
the pattern for all source positions. The phase determination of 
each of the spectral peaks allows resolving to finer and finer 
degrees as with a vernier caliper. 

(7) A convenient mathematical encryption of the spectral peak 
frequency and phase data (without using the final interpolation 
routine) uses a 4 digit number (base 4), where each digit defines 
one of four 90° positions for that spectral component. The MSB 
(Nf 0 ) corresponds to the observed fundamental frequency and the LSB 
(Nf 2 ) corresponds to the highest frequency component. Using D to 
represent the detector length, the pattern position (Y d ) is 
calculated as: 
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Y d = [(N f0 *4 3 + N fl *4 2 + N f2 *4 1 + N f3 *4°) / (4*f 3 ) ] *D 
The interpolated residual of the f 4 phase is included in this 
result to give a more accurate position along the detector surface. 
The pattern magnification is then determined by dividing the 
maximum detected single frequency (F d equals f 3 in this example) by 
the known maximum mask frequency (F m ) . 
= F d /F m 

The center of the mask segment illustrated by the point source (Y m ) 
is computed from the detector position and the pattern 
magnification by: 

The angle of incidence (0i) is determined from the mask height X m 
and the mask position by: 
9i = TaiT 1 (Y tt /XJ . 

(8) In the second variation, a variable frequency component is 
employed in addition to the multiple single frequency components. 
As seen in Figure 16, after magnification of the point source light 
beam image on the detector has been computed, the variable 
frequency component peak is calculated to yield coarse position 
data. Thereafter, in conformance with the first variation, phases 
of frequency components are determined, and so on. 

(9) Another aspect of the present invention is to determine 
the distance of the point source from the mask surface. This can 
be determined in two degrees of position. 

X s = X/d^-l) 



Y s = Y m 4- X s *Tan(9i) . 

The accuracy of the distance calculation depends on distance 
and detector size. In general, the accuracy will only be good at 
close distances, e.g., K m > 1.1, but for certain applications, this 
5 may be adequate. A second detector of this type placed in an 
orthogonal position to the first detector allows determination of 
a second radiation angle, as well as a distance projected 
orthogonally to its axis. From this information, the point source 
Q position (x,y, z) can be determined. 

%Q (10) In the second embodiment of the present invention, a 

% composite image is also used except that only one single frequency 
component is used in addition to a non-periodic function provided 
with unique auto-correlation. The single frequency function is 
used to determine the aperture image magnification. This scale 
factor is then used to resize the aperture image to correspond to 
the received image. 

(11) The two functions can then be cross-correlated in a 
manner explained in greater detail below to yield the precise image 
displacement on the detector surface. Generally speaking, only the 
20 function of the masked pattern that has a unique auto-correlation 
is employed in the correlating function. However, the single 
frequency component may also be used in the correlating function 
since this results in a periodic ripple superimposed on the unique 
cross-correlation function. In this way, the single frequency 
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component is only required for scale determination and is not 
required in the cross-correlation function. 

As such, it is a first object of the present invention to 
provide embodiments of a range adaptable system for determining 
angular position and distance of a radiating point source. 

It is a further object of the present invention to provide 
such a system in which a patterned mask or coded aperture is 
located at a distance from a linear array detector and interposed 
between the linear array detector and a point source of radiation. 

It is a still further object of the present invention to 
provide such a system which permits actual measurements of angular 
position and distance of a radiating point source from a detector 

Q without resort to comparison with reference data. 

m 

U It is a yet further object of the present invention to provide 

§5 such a system and method in which multiple single frequency 
components and phase methodology are employed in the measuring 
process . 

It is a yet further object of the present invention to provide 
such a system and method in which multiple single frequency 
20 components as well as a variable frequency component are employed 
in the measuring process. 

It is a still further object of the present invention to 
provide such a system and method, in another embodiment, in which 
only a single frequency component is used in addition to a non- 
25 periodic function having unique auto-correlation. 
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These and other objects, aspects and features of the present 
invention will be better understood from the following detailed 
description of the preferred embodiments when read in conjunction 
with the appended drawing figures. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 shows a schematic perspective representation of the 
spatial relationship between a point source of radiation, a window 
mask pattern, and a detector. 

Figure 2 shows a front view of the Figure 1 representation. 

Figure 3 shows a side view of the Figure 1 representation. 

Figure 4 shows a top view of the Figure 1 representation. 

Figure 5 shows a schematic representation of an analog mask 
pattern for a first embodiment of the present invention. 

Figure 6 shows a schematic representation of a digital 
detector pattern for the first embodiment of the present invention. 

Figure 7 shows an enlarged view of superimposed analog and 
digital masks of the first embodiment. 

Figure 8 shows a fast Fourier transform of the analog pattern 
of the first embodiment. 

Figure 9 shows a fast Fourier transform of the digital pattern 
of the first embodiment. 

Figure 10 shows the interpolated fast Fourier transform of the 
analog pattern of the first embodiment. 
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Figure 11 shows the interpolated fast Fourier transform of a 
particular frequency component of the analog pattern of the first 
embodiment . 

Figure 12 shows a fast Fourier transform of the mask with a 
first chirped frequency. 

Figure 13 shows a fast Fourier transform of the mask with a 
second chirped frequency. 

Figure 14 shows a fast Fourier transform of the mask with a 
third chirped frequency. 

Figure 15 shows the mask used to produce the fast Fourier 
transform plots illustrated in Figures 12, 13 and 14. 

Figure 16 shows a flowchart of operation of the first 
embodiment . 

Figure 17 shows a schematic representation of a function 
jij> containing a composite frequency used to produce an analog or grey 
:W scale mask in a second embodiment. 

Figure 18 shows a digital representation of the graph of 
Figure 17. 

Figure 19 shows a schematic representation of a fast Fourier 
20 transform of the mask pattern illustrated in Figure 17. 

Figure 20 shows a schematic representation of a fast Fourier 
transform of the schematic representation of the mask illustrated 
in Figure 18. 
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Figure 21 shows the correlation function produced by 
correlating the image with the random portion of the mask 
illustrated in Figure 21. 

Figure 22 shows the correlation function produced by 
correlating the image with the random portion of the mask 
illustrated in Figure 22. 

Figure 23 shows an analog transmissivity mask used in the 
second embodiment where the single frequency component and pseudo- 
random component are sequential. 

Figure 24 shows a digital transmissivity mask used in the 
second embodiment where the single frequency component and pseudo- 
random component are sequential. 

Figure 25 shows a flowchart of operation of the second 
embodiment . 

SPECIFIC DESCRIPTION OF THE PREFERRED EMBODIMENTS 
A first embodiment of the present invention is disclosed in 
detail with reference to Figures 1-16. A primary aspect of the 
first embodiment of the present invention describes a method for 
determining the angular position of a point source of radiation 
with respect to a detector by examining the scales and shifts of 
periodic components of the projected image of a variable 
transmissivity mask. The preferred embodiment encodes the 
transmissivity mask with several frequencies that have logarithmic 
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separations and computes the frequency spectrum of the received 
detector pattern. 

With reference to Figures 1-4, the configuration of the system 
used to practice the method is seen. With reference, first, to 
Figure 1, the system is generally designated by the reference 
numeral 10 and includes a point source emitter 11 designed to emit 
a point source of light having a beam 13 that may travel toward a 
coded planar mask 15. The mask is shown as including a plurality 
of lines 17 thereon defining a transmissivity pattern varying in 
one dimension (in the Y-axis) . The beam 13 travels through the 
mask 15 and impinges on a planar detector array 19. The mask and 
detector preferably lie in parallel planes. Figure 1 shows the 
image of the mask 15 projected onto the detector array 19. 
Reference numeral 21 refers to a mask segment that is projected 

p| onto the detector array 19 as shown by reference numeral 23, with 

ftj 

the length of the segment designated by the reference numeral 23 
being larger than the length of the mask segment 21 to a degree 
based upon the distance from the emitter 11 to the detector array 
19. The distance H between the coded mask 15 and the detector 

20 array 19 has a bearing on the ratio between the length of the mask 
segment 21 and the length of the array segment 23. As such, when 
the system is being set up to practice the inventive method, the 
distance H is programmed into the computer (not shown) . The length 
of the mask segment 21 is shown by the letter "d" and the length of 

25 the array segment 23 is shown by the letter "D." 
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Figures 2, 3 and 4 show front, side and top views, 
respectively, of the system 10. 

Figure 5 shows a function containing the composite frequency 
function used to produce an analog or grey scale mask, while Figure 
6 shows the digital version of the same function. Figure 7 is an 
enlargement showing both analog and digital representations. 

Figure 8 is a Fourier Frequency Transform (FFT) of the mask 
pattern of Figure 5, showing the spectral components, while Figure 
9 is the FFT of the mask of Figure 6. The analog mask produces no 
harmonic content, which shows up as noise in frequency domain and 
is preferred over the digital mask. 

Figure 10 shows the interpolated FFT plot of the mask pattern 

p of Figure 5, showing in high detail the four spectral peaks. 

pi j 

y| Figure 11 is an enlarged view of the FFT plot of Figure 10 showing 

Ui 

H the magnitude and phase of a single spectral component. 

am? 

SI 

From FFT plots of Figures 8-11, the spectral peaks may be 
determined, which indicate the received pattern frequencies. The 
corresponding phases indicate the position of each component. 
Because the point source is not located at infinity, the detected 

20 image pattern is never equal to the mask image, but is always 
magnified, with magnifications as much as 1.05:1 to 1.1:1. The 
spectral peaks are determined to a finer resolution than that of a 
detector pixel using interpolation. A typical application might 
use a 4,096 element CCD array with a mask pattern containing four 

25 frequencies, including the fundamental frequency f 0 with period of 
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4, 096 pixels, f x = 4*f G , f 2 = 16*4, and f 3 = 64*f Q . The mask length 
(-8192 pixels) is longer than the detector length to ensure that 
the detector is fully covered by the pattern for all source 
positions. The phase determination of each of the spectral peaks, 
successively, allows resolving the solution to finer and finer 
degrees as with a vernier caliper. A convenient mathematical 
encryption of the spectral peak frequency and phase data (without 
using the final interpolation routine) uses a 4 digit number (base 
4), where each digit N fx defines one of four 90° positions for that 
spectral component. The MSB corresponds to the observed 

fundamental frequency and the LSB corresponds to the highest 

yd 

jjj' frequency component. The pattern position is calculated as 

U Y d = [(N f0 *4 3 + N fl *4 2 + N f2 *4 1 + N f3 *4°) / (4*f 3 ) ] * D 

The interpolated residual of the f 3 phase is added to this 

31 result to give a more accurate angular position. Figure 11 shows 
the interpolated Fast Fourier Transform (FFT) of a particular 
frequency component. The phase in the vicinity of the peak 
response is linear and the phase and frequency corresponding to the 
FFT magnitude peak yields the desired pattern shift. 

20 A variation in the first embodiment uses an additional 

component to the mask pattern, which is a sinusoid with spatially 
varying frequency. A preferred spatial dependency is linear 
dependence. The FFT of the image received by the detector array 
allows determination of the centroid of the varying frequency 

25 component and from this information, the approximate pattern 
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position and corresponding angle of incidence are known. From this 
approximate solution, the phases of the single frequency components 
can be determined in sequence as explained above* The advantage 
and novelty of this variation is that the mask can be made much 
larger than the fundamental period of one array length. This also 
allows for the mask to be kept at a much greater distance from the 
detector surface, increasing the angular sensitivity and the 
ability to manufacture the device. 

The following equation expresses the mask transmissivity 
function: 



Mask, :- 0.5 + 



. MOO „ \ f i-20 „ 1 f i-4 ^ , 

cos 2-n + cos 2-n + cos 2-xc + cos 
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% Figures 12, 13 and 14 are FFT plots showing shifts in the 

y variable frequency component for two spatial shifts of the received 
jj] mask. Figure 15 is the mask used to produce FFT plots of Figures 
13 and 14. 

15 To account for magnification of the mask image at close 

ranges, the received variable frequency component is scaled by the 
reciprocal of the magnification as determined from detection of 
single frequency spectral components. A prior knowledge of the 
mask frequency components is assumed. 

20 Another aspect of this invention determines the distance of 

the point source from the mask surface. This can be determined in 
two degrees of position. 
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X s = H/d^-l) 

Y s = Y m + X s * Tan(9i) 

The accuracy of the distance calculation depends on distance 
and detector size. In general, the accuracy will only be good at 
close distances, e.g., > 1.1, but for certain applications, this 
may be adequate. 

A second detector of this type placed in an orthogonal 
position to the first detector allows determination of a second 
radiation angle, as well as a distance projected orthogonally to 
its axis. From this information, the point source position (x,y, z) 
can be determined. 

General applications requiring longer range operation will use 
three detectors, at least two of which are substantially separated 
to define three intersecting source planes that define the source 
position. 

In the method of measuring the position of the point source 
emitter 11 with respect to the detector array 19, reference is made 
to Figure 16. As a first step, the image of the mask projected 
onto the detector array 19 is first recorded. With the emitter 11 
activated, the Fast Fourier Transform (FFT) of the image on the 
detector array 19 is computed. Next, the highest frequency F d is 
identified. Then, the degree of magnification of the segment 23 as 
compared to the segment 21 is calculated. Next, if the process 
will involve use of a variable frequency mask, the variable 
frequency peak is determined to yield the "coarse" position. If 
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the variable frequency mask is not being employed, the phase at F 0 
is determined to yield coarse position. In either case, once the 
coarse position has been determined, the phases of the frequency 
components F a to F n are determined, followed by determination of the 
5 total pattern shift as a sum of the individual frequency component 
phase shifts. Thereafter, the pattern shift of the magnified image 
Y d is determined and the actual unmagnified mask pattern shift is 
calculated. From the information determined and calculated, the 

LI angle of incidence 0 V is computed. 

w 

;0 A second embodiment of the present invention describes a 

method for determining the angular position of a point source of 
radiation with respect to a detector by examining the scales and 
shifts of periodic components of the projected image of a variable 
transmissivity mask. This embodiment uses the system depicted in 
S{5 Figures 1-4 and encodes the transmissivity mask with a single 
^ frequency component and a broadband component and computes the 

frequency spectrum of the received detector pattern. The single 
frequency component (fj is generally chosen to be near the maximum 
frequency recordable by the imaging device and/or the maximum 
20 frequency within the diffraction limit, so that the best 
magnification resolution may be obtained. The broadband function 
used for correlation (f corr ) may be represented by a closed form 
function such as a polynomial or trigonometric function or it may 
be a numerical function, such as a pseudo-random sequence or 
25 universal redundant array. Ideally, the function will have a 
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nearly perfect impulse in the correlation with near zero 
correlation away from the peak. 

The mask transmissivity function is represented as 

mask(nT) = sin(2nf m nT) + f corr (nT) 
The magnified image function may be expressed as 
image (nT) = sin (2nf m nT/k m ) + f corr (nT/kJ 
where f m is the mask frequency and fJK is the detector frequency. 

If the correlating function is a closed form expression, then 
the image function is computed directly at n/k^ intervals. If the 
correlating function is a numerical or pseudo-random sequence, 
interpolation is used to compute the image mask at n/k m intervals. 

Figure 17 shows a function containing the composite frequency 
function used to produce an analog or grey scale mask, while Figure 

18 shows the digital version of the same function. 

hi 

jSp Figure 19 is a FFT of the mask pattern of Figure 17, showing 

the spectral components, while Figure 20 is the FFT of the mask of 
Figure 18. The analog mask produces no harmonic content, which 
shows up as noise in frequency domain and is preferred over the 
digital mask. Ideally, the correlating function spectrum should 

20 not overlap the single frequency spectrum to allow the best 
estimate of the single frequency peak. The example function shown 
(pseudo-random sequence) does overlap, but the effect on the 
spectral peak of the single frequency is quite small. 

From the FFT plot, the single frequency spectral peak is 

25 determined and the scale factor is computed by 
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Letting g(nT) represent the scaled mask correlation function, 

g(nT) = f corr (nT/k m ) 
and h(nT) represent the detected image, the cross correlation 
(R g , h (i)) is most efficiently computed using the FFT method and is 
expressed as 

Rg #h (T) = FFr x {FFT{g(nT) } . FFT*{h(nT)}} 
where the notation FFT' 1 represents the inverse FFT and FFT* 
represents the complex conjugate of the FFT. The peak of the 
correlation function indicates the lag or lateral distance that the 
image has projected along the detector surface- Note that the 
periodic function (a sine function in this case) may be kept in the 
correlating function, but that this adds a periodic component to 
the correlation function. 

Figure 21 shows the cross-correlation of the composite analog 
mask function in Figure 17 with the random component of the analog 
mask function. Figure 22 shows cross-correlation of the composite 
digital mask function in Figure 18 with the random component of the 
digital mask function. The correlation function peaks are clearly 
visible, with the digital mask function producing more artifacts 
due to digitizing. 

An alternative variation of this embodiment uses the unique 
correlation function peak only to determine a coarse position, 
while the fine position is determined by the phase of the single 
frequency component as in the first embodiment. 
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Another variation to this embodiment encodes the single 
frequency component and pseudo-random component sequentially in 
transmissivity mask instead of superimposing them. Figures 23 and 
24 show a schematic representation of this variation. This 
variation allows the non-periodic function to occupy the central 
portion of the mask, while the periodic function occupies the 
remainder of the mask, generally extending beyond the detector 
width. In this way, the non-periodic function image is fully 
detected at all angles of incidence, while the periodic function is 
again for scale identification. The identification of pattern 
shift is performed using the correlation technique described above. 
The correlation peak may be used to determine pattern fine shift or 
m it may be used only to determine coarse pattern shift, while the 
U phase of the single frequency component is used to determine the 
% fine pattern. The transmissivity functions may again be either 

m 

analog or digital as previously shown. 

With reference to Figure 25, the practice of the method of the 
second embodiment will now be described. First, as before, the 
detector array image is recorded, the FFT of the image is computed, 

20 the highest frequency F d is identified, and the magnification D/d 
is computed. Thereafter, the method of operation of the second 
embodiment diverges from the method of the first embodiment. The 
mask pattern is magnified by a factor of K,,. Thereafter, the 
detector image is compared to the magnified mask pattern using the 

25 technique of cross-correlation. From this information, the pattern 
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shift of the magnified image Y d is determined . From the 
information previously obtained, the actual unmagnified mask 
pattern shift Y m is computed. Finally, the angle of incidence 0 Y 
is computed. 

As such, an invention has been disclosed in terms of preferred 
embodiments that fulfill each and every one of the objects of the 
invention as set forth above, and provide a new and useful method 
of employing a range adaptable system for determining the angular 
position and distance of a radiating point source, of great novelty 
and utility. 

W Of course, various changes, modifications and alterations in 

Mi 

WJ the teachings of the present invention may be contemplated by those 

If! 

* skilled in the art without departing from the intended spirit and 

scope thereof. 

As such, it is intended that the present invention only be 
limited by the terms of the appended claims. 
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